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> . 1. b quarks in A-\-A collisions 

in 

In relativistic nucleus- nucleus collisions, heavy quarks (c,b) are produced in the early 
stage [|2], at time scales £~0.2-0.5 fm/c, in comparison to typical pion freeze-out 
time scales of i~10 fm/c. In particular, b quarks are interesting for two reasons: 
■ (1) The mass of the c quark m(c)=1.40 GeV is about the same order of magnitude 

as Aqcd— 1 GeV, but m(6)=4.75 GeV is significantly larger. Thus, for the b quarks, 
pQCD based cross section predictions are more reliable than for c quarks. 
\ (2) For Au+Au collisions at RHIC, the primordial temperature might be high 

enough to create cc pairs thermally, but not bb pairs. For t=0.6 fm/c, a temperature 
of T~0.3 GeV can be estimated At such a high temperature, QCD predicts that 
gluons, which are massless at T=0, acquire an effective mass, which turns out to be 
in the same order of magnitude as the temperature itself, i.e. m g ; uon ~0.3 GeV [0]. 
Thus the production process gg—>cc might be enhanced thermally. However, for bb 
production, the Boltzmann suppression factor is <bb> / <cc>=exp[— 2m(6)/2m(c)] 
; i ■ ~l/23, thus any thermal bb enhancement is unlikely. Thus, b quarks can be regarded 

as a pure probe of the early stage of the collision. However, bb events are rare. At 
RHIC, bb measurements are now feasible for first time, due to three reasons: (1) 
The bb production cross section at y/s=200 GeV (RHIC) is about a factor of ^10 3 
higher than at ^=18 GeV (SPS). (2) High collider luminosity provides a recorded 
data sample of 10 6 -10 7 events. (3) Recent technological and price development of 
fast CPUs created the basis for the possibility of event reconstruction in realtime 
for application as high level trigger on rare events. 

This paper describes the first attempt to observe bb production with the STAR 
experiment at RHIC. The main STAR detector is a large scale (R—2 m, L=4 m) 
TPC (Time Projection Chamber) [Q] with 5xl0 7 position pixels, each measuring 
the particle ionization dE/dx using a 10-bit ADC. For triggering on rare bb events, 
a Level-3 trigger CPU farm was used, described in Section 0] 

The 66-bound state T is specifically interesting, as the formation and re-breakup 
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rates are sensitive to color screening and thus the initial gluon density. STAR is 
particularly suited for the search for T production due to its large midrapidity 
coverage. Because of its high mass, the T is produced almost at rest, disfavouring 
gluon-boosted forward/backward emission of the decay products. 



2. The case of the T 



In the following, the term T will be used to represent the three S states Y(IS'), 
T(25) and T(35). T properties are listed in Tab. 1. The production ratio is 
independent from y/s in leading order QCD, thus ratios measured at y/s—\.% TeV 
are applicable to RHIC energies (-^5=200 GeV). 

Table 1. Mass m production ratio 
for midrapidity R pro d [El an d branching 
fraction [ |S] into e + e~, for the lowest 
three T states. 





m/GcV 


Rprod 


BR(e+e-) 


T(15) 
T(2S) 
T(3S) 


9.460 
10.023 
10.355 


72% 
18% 
10% 


2.38% 
1.18% 
1.81% 



The T(4s) and T(5s) shall not be considered here, as they are above the BB 
and B S B S thresholds, respectively. All T singlet 1 S, 1 P states are still unobserved, 
although they are expected to exist. So far, only T triplet states have been observed, 
because most data were recorded from e + e~ annihilation ( J PC ("f*) = l~~). Recently 
the discovery of T(1D) was reported [ [7], the only long-lived L=2 meson ever 
observed. The T decays into e + e~ with a branching fraction of ~2% (Tab. 1). 
Experimentally, this is the easiest exclusive final state to measure, whereas the 
decay into a ggg final state (3 jets) with the largest branching fraction of 77.4% is 
experimentally unobservable. At y / s=200 GeV, ~90% of the T are produced by 
gg^bb, ~10% by qq-^bb (q denoting u, d, s or c). As the total cross section for bb 
production in p+p is a(bb)=1.5 ub, one expects about one bb pair in ~52 Au+Au 
central events. With a production cross section of cr(T)=8.4 nb, it follows, that 
1/178 of bb pairs form a resonance 1 . A pure statistical production model, which 
works well for lighter particles, does not seem to be applicable For a rate 

estimate, the baseline prediction for cr pp (T) is taken from [ 9 . The cross section 
has a typical threshold behaviour a—a (l~ mj-/y/s) n , rising steeply in the RHIC 
energy region. Feed-down contributions are included in the calculation, e.g. only 
52% of the T(15) are produced directly. Feed-down contributions are 10% from 
the T(2S), 2% from the T(35), 26% from the XbO-P), and 10% from the x b (2P). 
With an absorption correction a, the fraction of the geometrical cross section f geo 
(10% for central collisions) and the fraction of hard processes fhard, the number of 
expected T per one central Au+Au collision at -^5=200 GeV can be estimated as 



N Au+Au 



Ay ■ s rec ■ (AA) a ■ f, 



hard 



J geo ' fgeo 



Bi?(e+O-^fl,=0 

2- 14.5%. (197.197)^.0.4 
F 7.26-0.1 



(1) 



As a comparison, 1/116 cc pairs form a J/ij). 
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3. Lepton Identification 

Electrons were identified by their specific ionization dE/dx in the TPC, which sat- 
urates at dE/dx~17-18 keV/cm for p>l GeV/c. Due to the mass dependance 
in the relativistic Bethe-Bloch formula, charged pions show increasing dE/dx~ll- 
14 keV/cm for p=l-10 GeV/c. As shown in Fig. [3 for this analysis electron candi- 
dates were identified by a cut on dE/dx: M1 , 




Fig. 1. TPC ionization dfi^ctx *1or cl Level-3 trigger (left) and offline (right) for 
Pt>3 GeV. For electron identification, a cut dE/dx>15 keV/cm was confirmed. 

The observed e/ir ratio for events with two p>3 GeV/c tracks is ~l/3.5, deter- 
mined from a two-gaussian fit (line in Fig. 1). The pion contamination in the 
electron sample is 17%. Imperfection of the fit in the region dE/dx~14-16 keV/cm 
is resulting from muon contamination, which is not taken into account for this pa- 
per 2 . Additional contamination from other particles (e.g. deuterons) is included 
in the systematic error. The electron identification was additionally checked with 
a partial electromagnetic calorimeter [ ^] (400 towers, Aip—60°, 0<rj<l). For 
dE/dx-selected electron candidates p/E=1.0±0.2 was achieved. 



4. The Level-3 Trigger 

The STAR Level-3 trigger performs full-event reconstruction in realtime, utilizing 
both fast CPUs and a fast network. 432 Intel i960 CPUs perform TPC cluster 
finding, 48 ALPHA DS-10 CPUs perform TPC track finding, and 3 Intel Pen- 
tium CPUs apply global algorithms, such as invariant mass calculations. For the 
10% most central Au+Au events (~130.000 clusters, ~4500 tracks), a sustained 
Level-3 processing rate of i?<40 Hz was achieved. For a solenoidal magnetic field of 
5=0.25 T, the pr resolution for (electron) tracks with pr=5 GeV/c was determined 
offline to be Apr/pr— 9% and ~17% (with and without vertex constraint, respec- 
tively), and, on the Level-3 trigger, Ap T /p T ~10% and -26% [IT. For B=0.5 T, 
all resolutions are to be multiplied by a factor 1/3. The dE/dx resolution is offline 
~8%, and on the Level-3 trigger ~11%. Further details are described elsewhere [ 

2 For dilepton invariant mass in the high mass region, mass differences between electrons and 
muons are negligible, and thus, for the muon contamination, an identical invariant mass slope is 
expected. 
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111) . For the T, the Level-3 trigger algorithm performed a trigger YES/NO decision, 
based upon an invariant mass cut for pt>3 GeV/c and dE/dx>15 keV/cm selected 
electron candidates. The T enhancement factor varied in the range 1.8-4.1, based 
on momentary accelerator luminosity conditions. 

5. Results 

The analyzed data sample is A r £t)e „ f =2,400,060 Au+Au 10% most central, addi- 
tionally iV B . ueTlt =212,030 Au+Au 10% most central, triggered by the Level-3 T 
algorithm, and A f E lJent =2,309,063 Au+Au minimum bias events. The data were 
taken during a period of ~100 days at y/s — 200 GeV. Assuming a number of 
1150 binary p+p collisions per one Au+Au collision, the data set corresponds to 
~4.2-10 9 p+p collisions. As a comparison, the amount of raw data on tape cor- 
respond to ~10xtotal data sample of a recorded 100 fb _1 B meson factory data 
sample. Fig.^shows the experimental results for the reconstructed invariant mass. 
A pseudorapidity cut of \rj\<l is applied. The top row shows di-electrons for unlikc- 
sign (left) and like-sign 3 (right) charge pairs. The straight line represents an ex- 
ponential fit to the m<8 GeV region, in order to visualize slope changes. The 
second row shows invariant mass for a control sample, i.e. 7r H 7r H candidates 
selected by dE/dx<15 keV/cm and normalized to the number of di-lepton candi- 
dates. The systematic error 4 consists of of a pr dependant Apt/pt resolution, a 
charged pion contamination of 17%, and an additional 10% systematic uncertainty 
induced by cut selection. The experimental data are compared to PYTHIA 6.212 
1121 calculations for the different subprocesses pp^rBB^e + vXe~vX (third row, 
left), pp^DD^e~isXe + VX (third row, right), Drell-Yan pp^e + e~ X (fourth row, 
left), and T production pp^>Y(lS,2S,3S)-^e + e~ (fourth row, right). The expected 
STAR mass resolution is Am=340 MeV for the T(15) alone, and Am=610 MeV 
for the T(1S,2S,35) combined (non-resolvable). The PYTHIA fc-factor 5 was set 
to fc=3.0. GRV94L was used as the parton distribution. All simulated mass dis- 
tributions are normalized to the experimental data sample (Section , taking into 
account enhancement factors on the Level-3 trigger. Two conclusions can be drawn. 
(1) The dilepton invariant mass distributions follow an exponential slope for m<9 
GeV, whereas the fit line is in the following refered to as reference line. An enhance- 
ment above the reference is visible for m>9 GeV, for both unlike-sign and like-sign 
di-electrons. The normalized di-pion control sample shows a non-exponential slope, 
for m>6 GeV beneath the reference line. As seen from reference line comparison, 
any possible contamination from the di-pions cannot explain the observed enhance- 
ment of the dileptons for m>9 GeV. From comparison to PYTHIA, it can be con- 
cluded that the major distribution to the di-electrons results from BB production 
with subsequent semileptonic decays. Drell-Yan docs not contribute significantly, 

3 The like-sign histogram is normalized according 2- V A r ++ ■ N . 
4 A11 error components were summed in quadrature. 

5 fc-factor is defined as ratio of cross sections for next-to- leading-order QCD 
and leading-order QCD. No higher loops were taken into account. 
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as the cross section is small a(qq— >e + e~)<10 nb. DD events are highly suppressed 
with a factor ~10 7 due to the electron p>3.0 GeV/c cut. 

(2) The unlike-sign dilepton distribution seems to indicate peak-like structure in the 
mass region m>9 GeV. In order to exclude possible accidentally used B=0 T events, 
which could fake high-p^ tracks, a second independent analysis confirmed the peak 
structure. The peaks remain visible when using a 2-dim (Level-3 trigger) or 3-dim 
(offline) track finder algorithm, or when perform track finding with or without vertex 
constraint. If the peak structure was induced by T production, it would indicate 
the unlikely fact, that the production cross section would be enhanced by a factor 
~10 3 compared to the baseline prediction (Section [3J. Any additional processes 
generating high mass di-leptons are not known. 



6. Systematic studies of a pp (T) cross section prediction 

The question of whether the observed peak structure might be induced by enhanced 
T production was investigated by systematic variation of three different basic QCD 
parameters in PYTHIA 6.212. 

(1) In the Cpp(T) calculation a fixed mass m(6)=4.75 GeV is assumed, but m(b) de- 
pends on the momentum transfer Q 2 in gg,qq^bb, i.e. m(6)~4.5 GeV for Q—l GeV/c, 
m(6)~3.5 GeV for Q=10 GeV/c, m(fe)~3.0 GeV for Q=100 GeV/c. 

(2) <7pp(T) also depends on the squared T wave function at the origin |4'(r=0)| 2 , 
which is calculated from the leptonic decay width 

r / T f+? ^_ 647r ag m e 2 |tt(r = 0)| 2 

e^e 

with exp. errors ±4%(ls),±16%(2s),±18%(3s). 

(3) a pp (T) depends on the strong coupling constant a s , which itself depends on 
the gluon distribution xG(x,Q 2 (x)). As an example, the early stage of an A+A 
collision might be governed by gluon saturation, simply given by the geometrical 
condition, that the probability of two colliding gluons in a gluon cascade is equal 
to one. 

3n 2 a s A x-G{x,Q 2 {x)) _ 

9T)2 X Td2 - 1 \ 6 ) 

^Saturation " ""A 

Additionally, a s is a function of temperature a s (T)~l/ln(T) [ 9,. At a scale of 
yG=V3GeV, [El quotes a s =0.280±0. 035 (expt)± 0.050 (syst)^;^ (theory). The 
total (quadratically summed) error is ±0.069, thus Aas/ag=25%. 
The variation of Aa pp (T) / a pp (T) from the baseline prediction (Ch.[5J| was investi- 
gated using PYTHIA 6.212, and found to be a factor 1.0—3.0 (only an increasing 
effect observed) for m(b), a factor 0.7—1.3 for |^I>(r=0)| 2 , and a surprisingly large 
factor of 10 —1.0 (only a decreasing effect observed) for as- However, in order 
to explain the observed peak structure in Fig. Q] (top, left) a <7 p)9 (T)-enhancement 
of a factor ~10 3 would be needed. Thus, at this stage, T production must be ex- 
cluded as explanation. For alternative explanations, two additional effects should 
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be considered: 

(1) A peak structure might be generated by BB (?ti=10.56 GeV) or B S B S (m=10.74 
GeV) threshold effects. Both meson species decay with large branching fractions 
of 10.5% and 8.1%, respectively. However, due to the undetected v in each decay 
and additionally hard-to-calculate QCD threshold dynamics, the size of the effect 
is difficult to estimate. 

(2) A B meson oscillates by second-order weak interaction into a B according to 
cos(Amt), with Am=0.503 ps" 1 [El- Using a B meson lifetime of ct=462 fim [|S] 
and a typical B meson Lorentz boost of 7=1.2 (corresponds to Eki n (B)=l GeV) it 
can be estimated for STAR, that ~2/3 of all B remain a B, and ~l/3 oscillate into 
a B. Every oscillation flips the sign of the lepton charge in the semileptonic decay, 
thus the oscillation could be observed in an unlike-sign/like-sign charge asymmetry 
vs. B meson decay vertex. Unfortunately, the STAR z vertex resolution has not 
been sufficient so far to resolve the oscillation. 

7. Apex 

Explanations of the observed enhancement and peak structures in the electron in- 
variant mass distributions might only be possible with a largely increased Au+Au 
data sample, which will be collected in 2004. Lepton identification at STAR will be 
improved by an electromagnetic calorimeter with |?y|<2 and Aip=2n coverage. For 
the envisaged RHIC-2 upgrade, starting possibly from 2007, a luminosity increase 
for Au+Au by a factor 40 to £=8xl0 27 cm _2 s -1 is planned, thus the T produc- 
tion will not be a rare process anymore. Additionally a new 2-layer pixel detector 
is planned at STAR, using a total of 9x 10 7 pixels of 20 fim 2 in area. Such a device 
will clearly enable an excellent measurement of B meson secondary vertices. 
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Fig. 1. Invariant mass for di-electrons and di-pions (control sample), for exper- 
imental data and PYTHIA 6.212. Statistical errors are shown as error bars, the 
white error band represents the systematic error. For details see Section |S] Please 
notice different y-axis scale by factor 10 -5 for T (bottom, right). 



